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ABSTRACT 

Context. Abell 1995 is a puzzling galaxy cluster hosting a powerful radio halo, but it has not yet been recognized as a obvious cluster 
merger, as usually expected for clusters with diffuse radio emission. 

Aims. We aim at an exhaustive analysis of the internal structure of Abell 1995 to verify if this cluster is really dynamically relaxed, 
as reported in previous studies. 

Methods. We base our analysis on new and archival spectroscopic and photometric data for 126 galaxies in the field of Abell 1995. 
The study of the hot intracluster medium was performed on X-ray archival data. 

Results. Based on 87 fiducial cluster members, we have computed the average cluster redshift (z) = 0.322 and the global radial 
velocity dispersion cr v ~ 1300 km s -1 . We detect two main optical subclusters separated by 1.5' that cause the known NE-SW 
elongation of the galaxy distribution and a significant velocity gradient in the same direction. As for the X-ray analysis, we confirm 
that the intracluster medium is mildly elongated, but we also detect three X-ray peaks. Two X-ray peaks are offset with respect to 
the two galaxy peaks and lie between them, thus suggesting a bimodal merger caught in a phase of post core-core passage. The third 
X-ray peak lies between the NE galaxy peak and a third, minor galaxy peak suggesting a more complex merger. The difficulty of 
separating the two main systems leads to a large uncertainty on the line-of-sight (LOS) velocity separation and the system mass: 
AV r f,Los = 600-2000 km s _1 and M sys = 2-5 xlO 15 M Q , respectively. Simple analytical arguments suggest a merging scenario for 
Abell 1995, where two main subsystems are seen just after the collision with an intermediate projection angle. 
Conclusions. The high mass of Abell 1995 and the evidence of merging suggest it is not atypical among clusters with known radio 
halos. Interestingly, our findings reinforce the previous evidence for the peculiar dichotomy between the dark matter and galaxy 
distributions observed in this cluster. 
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1. Introduction 

In the past decades, multiwavelength observations from ground 
and from space have dramatically shown the complexity of the 
physical phenomena occurring in galaxy clusters. An intrigu- 
ing aspect of these observations is the discovery of a growing 
number of clusters exhibiting diffuse radio emission (on Mpc 
scale), i.e. large-scale areas of radio emission without any ob- 
vious galaxy counterpart (Giovannini & Feretti 2002; Ferrari et 
al. 2008 ; Venturi 201 1 ). Particularly prominent are the radio fea- 
tures known as radio halos, which usually pervade the central 
cluster regions in a similar way to the intracluster medium (here- 
after ICM). Instead, radio emission areas found at the edges of 
clusters are known as radio relics. 
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The cause of radio halos and relics is still under investiga- 
tion. They are likely to result from synchrotron nonthermal radi- 
ation originating from relativistic electrons of the ICM moving 
in large-scale cluster magnetic fields. From a theoretical point 
of view, cluster mergers have been proposed as the key process 
for shedding light on the origin of these diffuse radio sources. In 
fact, the huge energy of these events could reaccelerate mildly 
relativistic electrons to relativistic energies and amplify the clus- 
ter magnetic fields (e.g., Feretti 1999). In particular, radio relics 
seem to be connected with large-scale shock waves occurring 
during mergers (e.g., Ensslin et al. 1998; Hoeft et al. [2004). 
Instead, radio halos are probably related to the turbulent mo- 
tions of the ICM following a merger (e.g., Cassano et al. 2006; 
Brunetti et al. 2009), but the precise scenario is still being de- 
bated. 

X-ray observations have been crucial to deriving the dy- 
namical state of clusters hosting diffuse radio emission. Several 
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statistical studies (see, e.g., Buote [20021 Cassano et al. 2010; 
Rossetti et al. 1201 11 1 have discovered interesting correlations be- 
tween the properties of radio halos and relics and the ICM X- 
ray luminosity and temperature (Giovannini & Feretti 2002] and 
refs. therein). This is also true when comparing point-to-point 
the X-ray and radio surface brightnesses (Govoni et al. 2001). 
Nevertheless, in a pilot study using the Sunyaev-Zel'dovich ef- 
fect, Basu (120121 1 finds the lack of bimodality in the radio power 
- integrated SZ effect measurement diagram, the contrary of 
what is found in the radio power - X-ray luminosity diagram 
(Brunetti et al. 2007). This study shows the need to adopt more 
investigation techniques in addition to the X-ray data analysis. 

Optical observations can be very helpful when checking for 
mergers in clusters with diffuse radio emission and studying their 
internal dynamics, too (e.g., Girardi & Biviano 2002). In partic- 
ular, combined X-ray/optical studies can be very effective at re- 
vealing and quantifying the level of substructure, when checking 
for premerging subsystems and/or merger remnants. The power 
of this approach comes from the fact that mergers affect the ICM 
and the galaxy distributions in different ways, as shown by nu- 
merical simulations (e.g., Roettiger et al. 1 1 9971 >. Thus, X-ray and 
optical observations complement each other to provide a more 
complete picture of merger events. 

It is with this scientific rationale in mind that we have be- 
gun a long-term optical observational program to investigate 
the properties of clusters hosting radio halos and/or relics: the 
DARC ("dynamical analysis of radio clusters") program (see 
Girardi et al. 2010). Among the dozens of clusters with known 
diffuse radio sources, we decided to perform an optical and X- 
ray investigation of the interesting cluster Abell 1995 (hereafter 
A1995). 

In the X-ray band, A 1995 appears as a luminous and hot 
cl uster: L x (0. 1-2.4 keV)=13.42xl0 44 /z 50 2 erg s" 1 (Bohringer et 
al. l2000l l; kT x = 7-9 ke V (from Chandra data, see e.g. J3aldi et 
al. 120071 Bonamente et al. l2008l and Ehlert & Ulmert 120091 . 

At optical wavelengths, A 1995 is a rich cluster (Abell rich- 
ness class = 1; Abell et al. 1989). Its light distribution is quite 
elongated in the NE-SW direction, but the mass shows a more 
circular distribution and is more concentrated than galaxies, as 
shown by the weak gravitational lensing reconstruction by Dahle 
et al. ( 120021 1. Holhjem et al. (2009) confirm this interesting dis- 
crepancy between light and mass distribution (see discussion in 
Sect. |5]l. 

As for previous redshift data, Patel et al. (2000) obtained 
spectra for 15 member galaxies and estimate a cluster redshift of 
Z = 0.322+0.001 and a radial velocity dispersion cr v = 1282+J20 

km s -1 . Irgens et al. (2002) confirm these measurements on the 
basis of 20 (unpublished) redshifts, six of them in common with 
those of Patel et al. (EUOOl . with z = 0.3207 + 0.0001 and 
cry = 1130^q km s _1 . Moreover, they find good agreement 
between the galaxy velocity dispersion and the dark matter ve- 
locity dispersion obtained from the weak gravitational lensing 
analysis (cr DM = 1240 + 80 km s~'). 

The X-ray ROSAT-HRI emission is peaked on a central 
bright galaxy and shows modest elongation in the NE-SW di- 
rection, which is not clearly separated from the emission of two 
very bright pointlike sources (see Fig. 1 of Patel et al. 2000, see 
also Ota & MitsudaHxHJl. Ota & Mitsuda d2004t classify A1995 
as a regular cluster from the stability of the X-ray centroid po- 
sition and the good fit with a single y6-model profile. However, 
despite its regular appearance, A 1995 has a very large cooling 
time (f coo i = 10.7 Gyr; see Fig. 12 of Ota & Mitsuda [20041 and 
thus no evidence of a cool core. 



Summarizing previous results, there is some hint but no clear 
evidence of substructure in A1995, and indeed, on the basis of 
optical and X-ray appearance, Pedersen & Dahle (2007 ) include 
this cluster in the sample of relaxed systems. 

About the radio wavelengths, Owen et al. d 19991 1 first re- 
ported a possible detection of a diffuse radio source in A 1995. 
Giovannini et al. (2009 ) analyzed new VLA data and discovered 
an evident radio halo in this cluster, with a size of ~ 0.8 fci 
Mpc and radio power of P l4 ghz = 1.3 X 10 24 W Hz -1 . The 
radio halo appears somewhat elongated in the NE-SW direction 
(Giovannini et al.[2p09i see also Fig. [I). 

In the context of our DARC program, we proposed new spec- 
troscopic observations of A 1995 with the Telescopio Nazionale 
Galileo (TNG). We also performed new photometric observa- 
tions at the Isaac Newton Telescope (INT) and used archival 
data of the Sloan Digital Sky Survey (SDSS). As for the anal- 
ysis in the X-ray band, we used archival data downloaded from 
the Chandra Archive. 

In this paper, Sect. 2 presents the new spectroscopic and pho- 
tometric data. Section 3 describes the analysis of the optical data, 
while Sect. 4 presents the analysis of the Chandra archival data. 
Finally, in Sect. 5, we discuss our results and propose a scenario 
for the dynamical status of A 1995. 

Throughout this paper, we adopt a flat cosmology with Ho = 
70 km s _1 Mpc" 1 (h 70 = H /1Q km s" 1 Mpc" 1 ), Q = 0.3 and 
Qa = 0.7. With this choice of the cosmological parameters, the 
scale is ~ 280 hzl kpc/arcmin at the redshift of A1995. Errors 
are given at the 68% confidence level (hereafter c.l.), unless oth- 
erwise stated. 



2. Galaxy data and catalog 

2.1. Spectroscopic observations 

We performed spectroscopic observations of A 1995 in May 
2009. As usual for the clusters in our DARC sample, we used 
the instrument DOLORES @TNC3 in MOS mode with theLR- 
B grism, which covers the wavelength range 3000-8430 A. In 
summary, we obtained 143 spectra from four observed masks. 
For each mask, the total exposure time was 5400 s. 

Reduction of spectra and radial velocities computation with 
the cross-correlation technique (CC; Tonry & Davis [19791 were 
performed using standard IRAF3 tasks, as for other clusters in- 
cluded in our DARC sample (for a detailed description see, e.g., 
Boschin et al. 2012 , hereafter B12). In eight cases (IDs. 02, 03, 
08, 14, 16, 29, 52, and 91; see TablelB the redshift was estimated 
with the EMSAO package (based on the wavelength location of 
emission lines in the spectra). Our spectroscopic catalog lists a 
total of 126 galaxies. 

After comparing the velocity measurements for galaxies ob- 
served with multiple masks (see discussion in, e.g., Boschin et 
al. 120041 Girardi et al. 1201 11 1, we corrected the velocity errors 
provided by the CC technique by multiplying them for a factor 
~2. After taking the above correction into account, the median 
value of the cz errors is 74 km s _1 . 

We also compared our data with those of Patel et al. (2000), 
finding 13 out of 15 galaxies in common. The two redshift mea- 
surements agree with a one-to-one relation, but the x 2 of the fit 
reaches a reasonable value only by multiplying their errors by at 
least a factor 1.5. This leads to their typical errors being three 



1 see http://www.tng.iac.es/instruments/lrs 

2 see http://iraf.net 
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Fig. 1. Optical/X-ray/radio view of A1995 (direct sky view with north up). A smoothed Chandra 0.3-7 keV image (orange and 
yellow colors) of A 1995 is superimposed on an /?H-band image taken with the INT. Thin contours are the radio contour levels from 
Giovannini et al. (2009; VLA data at 1.4 GHz with discrete sources subtracted, courtesy of F. Govoni). Thick contours are the mass 
distribution contours from Holhjem et al. (2009). Numbers highlight notable galaxies mentioned in the text. 



times larger than our typical errors. Considering these large er- 
rors, we preferred not to add this additional information and to 
study our (homogeneous) sample. 

2.2. Photometric observations 

We performed photometric observations of A 1995 in January 
2008 by using the Wide Field Camera (WFC0) mounted on the 
2.5m INT Telescope. The sky conditions were photometric and 



see http://www.ing.iac.es/Astronomy/instruments/wfc/ 



the seeing was ~1.4". In particular, we observed with the Harris 
Bh (15 exposures of 600 s) and Ru (18 exposures of 300 s) filters. 
This means a total exposure time of 9000 s and 5400 s in each 
band. 

We reduced the data and produced our photometric galaxy 
catalog by using standard procedures (see, e.g., B12 for details 
on the reduction of the WFC images). After transformation of 
Z?h and Rh magnitudes into the B and R Johnson-Cousins mag- 
nitudes (Johnson & Morgan 19531 Cousins 1 19761 1 and magni- 
tude correction for the galactic extinction (with Ab ~ 0.06 and 
Ar ~ 0.04, respectively; Schlegel et al. 1998), we estimated that 



3 



W. Boschin et al.: Structure of the cluster Abell 1995 



□ 



* 

. 



b , 



- t£> 




~* • 



. 



® 



00' 76 



* 



B. 



Q ^Jfft" 



« ® 





r 





* ; ^ v 



® 



• « 00 J- 
■ m » h 



1] 



□ e 



Fig. 2. Field of the cluster A1995 (direct sky view with north up) observed with the INT in the Ru-band. Cluster members are 
indicated by circles, while nonmember galaxies are shown by squares (see TableQ}. Numbers highlight the IDs of notable galaxies 
as in Fig.Q] 



the completeness of the photometric sample is 90% for R < 21.2 
and B < 22.7. 



which is a cluster member. Instead, 1452+5801 is likely a back- 
ground source. 



2.3. Galaxy catalog and notable galaxies 

Table Q] collects all the spectroscopic and photometric infor- 
mation for the 126 galaxies with a measured redshift (see also 
Fig. 0: ID and IDm (Cols. 1 and 2) are the identification num- 
ber of each galaxy and member galaxies, respectively; Col. 3 
reports the equatorial coordinates, a and 6 (J2000); Cols. 4 and 
5 list the B and R magnitudes; Col. 6 lists the radial (heliocentri- 
cally corrected) velocites, v = czq, with their errors, Av (Col. 7). 

A1995 hosts a dominant galaxy (ID. 61, R = 17.78, here- 
after BCG) ~ 0.8 mag brighter than the second brightest mem- 
ber galaxy (ID. 76). The X-ray ROSAT-HRI emission peaks on 
the BCG (e.g.. Patel et al. 120001 . 

There are several X-ray and radio emitting galaxies in the 
field of A1995. Among them, our ID. 89 (a background galaxy) 
is an evident pointlike source in Chandra archival data (see 
Fig. Q]). Cooray et al. (1998) highlight two bright radio point- 
like sources in A1995: 1453+5803 and 1452+5801 Takin g a 
look at the radio map provided by Giovannini et al. (2009] see 
their Fig. 8, right panel), we identify 1453+5803 with our ID. 71, 



3. Analysis of the optical data 

3.1. Member selection 

As for other DARC clusters, we selected cluster members by 
running two statistical tests. First, we used the 1D-DEDICA 
method (Pisani [T993I and [1996I I on the 126 galaxies with red- 
shifts. The method detects A1995 as a significant peak (at >99% 
c.l.) in the velocity distribution located at z ~ 0.322. The peak 
includes 94 (provisional) member galaxies (see Fig. [3}. 

Then, we used a second tool for member selection, which 
uses both the spatial and velocity information: the "shifting gap- 
per" method (Fadda et al. 1 19961 see also, e.g., Girardi et al. 1201 II 
for details on the application of this technique). Here, we only 
point out that the method needs the definition of a cluster cen- 
ter. For A1995, we chose the location of the BCG (see Table[TJ. 
The application of the "shifting gapper" rejected another seven 
galaxies in the outskirts of the cluster (Fig. [4]- top panel). Finally 
we obtained a sample of 87 fiducial cluster members (Fig. [5]). 
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Fig. 3. Histogram of 126 galaxy redshifts. The 94 (provisional) 
cluster members are highlighted with the solid line (see text). 
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Fig. 4. The 94 (provisional) cluster members (see also Fig. [3} 
in the phase space (see top panel). The ordinate is the rest- 
frame velocity, the abscissa the projected clustercentric distance. 
Galaxies rejected by the "shifting gapper" procedure are shown 
with (red) crosses. For the cluster center we adopt the location 
of the BCG (see text). Big green circles and small points show 
the differential and integral profiles of the mean velocity (in the 
middle panel) and the radial velocity dispersion (in the bottom 
panel). In the bottom panel, two horizontal lines mark the range 
of possible values for the ICM temperature (7-9 keV) with their 
respective errors transformed to cry (see Sect.[5]for details). 



Table 1. Catalog of 126 galaxies in the field of A1995 with mea- 
sured radial velocities, where f highlights the ID. of the BCG. 



ID 


IDm 


a,<5(J2000) 


B 


R 


V 

(km 


Av 


01 




14 5? 18 17 

1 T JZ lU.li., 


+58 01 33 3 

TJO V/1 JJ.J 


19.18 


18.16 


18112 


101 


02 




14 5? 74 17 

1^ J Z Z.T. 1 / , 


+58 00 77 5 
tjo z / . j 


22.31 


21.18 


74182 


82 


03 




14 5? 76 31 

1^ J Z Z\J. J 1 . 


+58 0? 34 

TJO V/Z JT.V/ 


20.50 


18.75 


70357 


170 


04 




14 5? SO R? 

1^ JZ JU.OZ., 


+58 0? 58 7 

TJO v/Z JO.Z 


19.45 


18.40 


49107 


129 


05 


1 


14 5? 33 0? 

1^ Jz JJ.UZ., 


+58 0? 77 6 

l JO V/Z Z / .V/ 


22.62 


20.09 


95649 


48 


06 


2 


14 5? 33 1 5 

1^ Jz JJ.l J, 


+58 00 1 3 9 

TJO V/ \J 1 J.7 


22.75 


20.38 


97643 


72 


07 




14 5? 33 19 

1^ Jz JJ.17, 


+57 59 33 7 

T J / J7 JJ.Z 


22.93 


20.35 


121114 


208 


08 


3 


14 5? S4 S9 


+57 59 54 6 

T J / J7 JT.VJ 


22.77 


21.41 


95562 


100 


09 




14 5? S4 Q? 


+58 00 40 5 

TJO v/V tv/. J 


22.54 


19.91 


128719 


141 


10 




14 5? 36 19 

it j z. jvj.i,?} 


+58 00 40 6 

TJO v/U tU.u 


21.77 


20.36 


103177 


105 


1 1 


4 


14 5? 36 53 

1^ J Z JU. J J, 


+58 00 44 9 

TJO v/U TT. y 


22.42 


20.11 


95879 


77 


12 




14 5? S7 80 

1^ JZ J / .OV/, 


+ 58 01 76 9 

TJO v/1 ZU.7 


21.99 


19.96 


138505 


73 


13 




14 5? 38 1 8 

1^ JZ JO.IO, 


+57 59 1 5 7 

T J / J7 1J.Z 


22.67 


20.78 


134283 


128 


14 




14 5? 39 1 ? 

1^ JZ J7.1Z, 


+58 00 54 

TJO v/U Jt.v/ 


22.82 


21.22 


143292 


100 


15 


5 


14 5? 39 87 


+58 0? 78 6 

TJO v/Z ZO.U 


23.42 


20.93 


94671 


73 


16 




14 5? 40 85 

1^ JZ tVJ.OJ, 


+57 59 19 3 

TJ/ J7 1 ~. J 


23.29 


21.42 


103596 


100 


17 




14 5? 40 89 

1^ JZ tU . O " ; 


+58 00 14 5 

TJO \J\J It. J 


22.36 


19.50 


128750 


72 


18 


6 


14 5? 41 04 

1^ JZ I 1 . V/T , 


+58 03 38 5 

TJO v/J J O.J 


22.07 


19.58 


95191 


41 


19 


7 


14 52 41.47 


+58 04 07 7 

TJO v/t VJZ. / 


22.62 


20.15 


96695 


97 


20 


8 


14 5? 4? 96 

It JZ Ti. . ZU; 


+57 59 79 3 

TJ/ J7 Z7.J 


23.54 


20.83 


96151 


76 


21 


9 


14 5? 4? 7? 

IT JZ Ti. . / Z., 


+58 00 41 7 

TJO V/V Tl.Z 


22.95 


20.43 


94533 


85 


22 




14 5? 4S S(S 

It Jz tJ.JU, 


+58 03 04 4 

TJO v/J V/t.t 


21.82 


19.85 


142517 


1 17 


23 


10 


14 5? 43 78 

IT JZ T J . / O, 


+58 01 79 5 

TJO v/1 Z7.J 


21.29 


18.81 


98437 


65 


24 


1 1 


14 5? 45 07 

It JZ tJ.W / . 


+58 03 1 3 3 

TJO v/J 1 J.J 


22.27 


19.72 


93027 


120 


25 


12 


14 5? 46 39 

It JZ tU.J", 


+58 01 17 4 

TJO V/1 1Z.T 


21.82 


19.50 


95659 


88 


26 


13 


14 5? 46 53 

IT JZ TVI.JJ, 


+58 01 36 7 

TJO v/1 JU. / 


23.31 


20.62 


97678 


116 


27 


14 


14 5? 47 96 

It JZ t / . v 1. 


+58 01 57 7 

TJO \J 1 J / .Z 


22.04 


19.45 


95027 


49 


28 


15 


14 59 47 45 

XT T/ .T^J* 


+58 01 37 

r V/1 J / .V/ 


22.52 


19.76 


98072 


108 


29 




14 S? 47 8fS 

It JZ t/.OU, 


01 10 

TJO Ul lu.u 


22.66 


21.46 


134784 


100 


30 


16 


14 52 47.86, 


+58 02 01.0 


21.73 


19.37 


93980 


45 


31 




14 52 48.86, 


+58 04 39.9 


20.17 


18.99 


20008 


68 


32 


17 


14 52 48.98, 


+58 01 48.2 


23.16 


20.54 


96159 


85 


33 


18 


14 52 49.20, 


+58 02 18.7 


22.50 


19.87 


95919 


85 


34 


19 


14 52 49.92, 


+58 00 26.3 


20.80 


19.18 


99092 


87 


35 


20 


14 52 50.11, 


+58 01 26.9 


21.37 


18.88 


94063 


57 


36 


21 


14 52 50.33, 


+58 01 50.7 


23.05 


20.32 


98429 


77 


37 




14 52 50.67, 


+58 04 06.1 


23.84 


20.76 


157123 


60 


38 


22 


14 52 51.17, 


+58 02 04.6 


22.37 


19.68 


97314 


81 


39 


23 


14 52 51.38, 


+58 02 36.5 


21.99 


19.84 


94876 


57 


40 


24 


14 52 51.94, 


+58 02 33.4 


22.00 


19.21 


95569 


49 


41 


25 


14 52 51.96, 


+58 02 10.1 


21.93 


19.29 


97901 


48 


42 




14 52 52.03, 


+58 00 45.9 


23.46 


21.36 


184420 


68 



3.2. Global cluster properties 

The first and second moments of a distribution can be efficiently 
computed by using the biweight estimators of location and scale 
(Beers et al. 1990). Their application to the velocity distribu- 
tion of our 87 cluster members provided a measurement of the 
mean cluster redshift and of the global radial velocity dispersion, 
where we found <z> = 0.3217+ 0.0005 (i.e. (v) = 96437±140 
km s~') and cry = 1302^j 7 km s _1 , respectively. 

Analysis of the velocity dispersion profile (Fig. HJ suggests 
that the value computed for cry is quite robust. In fact, the inte- 
gral profile is asymptotically flat within the errors. Instead, the 
decrease in the differential profile agrees with the fact that the 
cluster is free of interlopers, as also seen in the velocity vs. clus- 
tercentric distance diagram (Fig. [4]- top panel). As for the mean 
velocity profile, the modestly higher values of the mean velocity 
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Table 1. Continued. 



ID 


IDm 


a,<5(J2000) 


B 


R 


V 

(km 


Av 


43 


_ 


14 52 52.44, +58 00 49.1 


22.82 


21.37 


214244 


61 


44 


26 


14 52 52.54, +58 05 01.7 


22.33 


19.81 


97005 


61 


45 


27 


14 52 52.77, +58 04 39.0 


21.34 


19.95 


93754 


72 


46 


28 


14 52 52.87, +58 02 40.7 


20.89 


18.79 


96237 


40 


47 


29 


14 52 52.97, +58 03 00.8 


20.70 


18.81 


93425 


61 


48 


30 


14 52 53.31, +58 02 17.7 


21.96 


19.43 


93727 


49 


49 


31 


14 52 53.81, +57 59 33.9 


21.54 


19.24 


95494 


180 


50 


32 


14 52 53.81, +58 01 53.9 


20.28 


19.22 
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in the external region (Fig. [4]- middle panel) probably come from 
the slightly larger sampling in the NE external regions, where 
galaxies have higher velocities (see Sect. I3.51 >. 

3.3. Analysis of the velocity distribution 

Deviations of the velocity distribution from Gaussianity can pro- 
vide signs of disturbed dynamics. This can be checked by apply- 
ing classical shape estimators (Bird & Beers [19931 1. We did not 
not find significant evidence of departures from Gaussianity by 
using the skewness, the kurtosis, and the STI estimators. 

We also searched for significant gaps in the velocity distri- 
bution. In particular, we performed the weighted gap analysis 
(Beers et al. I1991l and [1992). We detected one significant gap 
that divides the cluster into two groups (see Fig . [5j : GV1 (with 
42 galaxies and lower velocities) and GV2 (with 45 galaxies and 



higher velocities). The BCG resides in the GV1 group, but it lies 
on the border with GV2. 

When considering the spatial distributions of the galaxies of 
GV1 and GV2, we found that they are different at the 99.94% 
c.l. according to the 2D KS-test (Fasano 1987, see our Fig. [6]). A 
statistical test useful to search for eventual subsets in the veloc- 
ity distribution is the lD-Kaye's mixture model test (1D-KMM; 
Ashman et al. 1994; see also, e.g., B12 for a description of the 
method). 

We applied this technique by assessing whether a two- 
Gaussian partition (accordingly to the detection of the two 
groups GV1 and GV2) can provide a significantly better fit to 
the velocity distribution than a sole Gaussian. The result of the 
1D-KMM test is negative. However, the best-fit result of the 
1D-KMM method gives two groups (KMM1D-1 and KMM1D- 
2) of 40 and 47 galaxies, very similar to the groups GV1 and 
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Fig. 5. The 87 galaxies recognized as cluster members. Top 
panel: radial velocity distribution with the arrow indicating the 
velocity of the BCG. Red and blue Gaussians are the best two- 
group partition fits according to the 1D-KMM test (see text). 
Bottom panel: Stripe density plot. The position of the significant 
gap is marked by an arrow. 
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Fig. 6. 2D distribution of 87 cluster members. Galaxies with 
smaller (larger) symbols have higher (lower) radial velocities. 
Blue circles and red crosses identify galaxies of GV1 and GV2. 
The origin of the coordinates is the location of the BCG. The 
solid and dashed lines indicate the position angle of the cluster 
velocity gradient (see text) and relative errors, respectively. 




4 2 -2 -4 

X [arcmin] 

Fig. 7. Isodensity contours of the cluster members spatial distri- 
bution computed with 2D-DEDICA. The origin of the coordi- 
nates is the location of the BCG (the big cross). Small circles 
mark the positions of cluster members relative to the BCG. 



GV2 (but note that the BCG is now belonging to KMM1D-2, 
the high-velocity group). The spatial distributions of the galaxies 
of KMM1D-1 and KMM1D-2 are different, too (at the 99.75% 
c.l.). In Fig. [5] and Table [2] we present the results for the two 
Gaussians with the velocity dispersions computed by the 1D- 
KMM procedure, which considers the membership probabilities 
of the galaxies to belong to a group, rather than the velocity dis- 
persions computed on the galaxy group populations after the as- 
signment. In this way we could bypass the artificial truncation 
of the velocity distribution tails of KMM1D-1 and KMM1D-2, 
thus minimizing the danger of heavily underestimating the ve- 
locity dispersions of the two subclusters. 

3.4. Analysis of the galaxy spatial distribution 

We applied the 2D-DEDICA method to the sky positions of 
A1995 member galaxies to search for eventual subsets in the 
galaxy spatial distribution. We found a NE-SW elongated struc- 
ture with two significant peaks, one of them (the NE one) cen- 
tered on the BCG (Fig. [7]). However, this finding could be af- 
fected by the magnitude incompleteness of our spectroscopic 
sample. To test the robustness of the 2D-DEDICA results, we 
resorted to the photometric catalogs. 

We used the color-magnitude relation (hereafter CMR) of 
early-type galaxies (the dominant galaxy population in the dens- 
est, internal cluster regions; e.g. Dressier 1980) to select likely 
cluster members from our photometric sample (see B12 for de- 
tails on the technique used for the determination of the CMR 
and the selection of member galaxies). We found B-R = 4.164 - 
0.081 x R (see Fig. [8). Figure [9] illustrates the contour map for 
the likely cluster members (513 with R <2 1 in the whole INT 
field): we confirm that A1995 is described well by an elongated 
structure. The two most significant peaks have similar densities, 
one centered on the BCG galaxy and one ~ 1.5'at SW. 

Table [3] lists information for these two main peaks (2D- 
NE and 2D-SW clumps), including Ns, the number of mem- 
bers (Col. 2), the peak densities, ps, relative to the densest peak 
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Table 2. Kinematical properties of several subsystems. 



System 



< V > 
km s _1 



km s -1 



M 



Notes 



Whole system 
KMM1D - 1 
KMM1D - 2 
KMM3D - 1 
KMM3D - 2 
HT1 
HT2 



87 96437 ± 140 
40 95130 3 
47 97472 a 
49 96055 ±200 
38 96850 ± 178 
51 95829+ 134 
16 98200 ± 90 



1302+ 1 , ? 
954 a 
805 a 

1080+}^ 
950 + ™ 
342 !2 



3.0 
1.2 
0.7 
3.6 
1.7 
1.2 
0.05 



ID low velocity subcluster 
ID high velocity subcluster 
3D low velocity subcluster 
3D high velocity subcluster 
main system 
secondary (high velocity) system 



Here we consider < v > and cry as given by the KMM software, where galaxies are opportunely weighted according to their membership 
probability (see Sect. l3.3l >. 
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Fig. 8. B-R vs. R diagram for the 126 galaxies with spectro- 
scopic information. Small black squares are cluster members, 
while big red squares represent field galaxies. The solid line 
shows the best-fit CMR as computed from cluster members; the 
dashed lines mark the region where likely cluster members were 
selected from the photometric catalogs. 




5 o 

X [arcmin] 

Fig. 9. Results of the 2D-DEDICA method (blue isodensity con- 
tour lines) applied to likely A 1995 members (see text) with 
R <21. The origin of the coordinates is the location of the BCG 
(the big cross). Small circles mark the positions of likely cluster 
members relative to the BCG. 



(Col. 4), and the significance of the peaks estimated through the 
value of x 2 (Col. 5). We also detect a minor peak (2D-NENE 
in Table [3] the third in galaxy density) whose statistical signif- 
icance, although nominally acceptable, is much lower than the 
two main peaks. 

We also considered the SDSS photometric catalogs (already 
corrected for Galactic absorption). In this case, we selected 
likely member galaxies by considering the CMRs in the (r'-i' 
vs. r') and (g'—r' vs. r') color-r' mag diagrams (see Goto et al. 
|2002]and B12). The two CMRs are r'-/'=1.009-0.022x r' and 
g'-r'=2.878-0.062x r' , respectively. The result for r' < 21.5 
(the limit for the SDSS star/galaxy classification) confirms the 
results reported above. 



3.5. 3D analysis: combining spatial and velocity information 

The 3D tests searching for correlations between positions and 
velocities of member galaxies are powerful tools for revealing 
real substructures in clusters. First, we checked for a velocity 
gradient in the set of the 87 fiducial cluster members (see, e.g., 
den Hartog & Katgert [19961 and Girardi et al. 119961) . We found 
a significant (at the 96% c.l.) velocity gradient with PA = 48^g 
degrees (counter-clock-wise from north). This means that the 
NE region of the cluster is populated by higher velocity galaxies 
(see Fig. [6j>. 

A classical test to detect the presence of substructures is the 
A-statistics (Dressier & Schectman l 19881 hereafter DS-test). We 
applied this test by defining, for each galaxy, the 5 parameter 
(determined considering the 7Y nn =10 neighbors of each galaxy; 
see, e.g., B12 for a description of the method) and computing 
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Table 3. Substructure from analysis of the INT photometric sam- 
ple. 



Subclump 



N & a, <5(J2000) p s X \ 
h : m : s ° : ' : " 



2D-SW(INTfi<21) 21 
2D - NE (INT R < 21) 23 
2D - NENE (INT R < 21) 13 



14 52 48.9 + 58 01 48 1.00 20 
14 52 57.6 + 58 03 03 0.99 18 
14 53 04.1 +58 03 50 0.70 11 



the departure A of the local kinematic quantities from the global 
parameters. We also applied the e-test (Bird 1993 ) and the a-test 
(West & Bothun 1990; see also, e.g., Pinkney et al. 119961 and 
Ferrari et al. 120031 for details of these 3D tests). Both the DS- 
test and the a-test detected the presence of substructures (at the 
98.5% and 95% al., respectively). 

A better interpretation of the results of the DS-test can be 
reached by splitting the contributions of the local mean velocity 
(estimator 6y) and dispersion (estimator <5 S ) to the classical 5 pa- 
rameter (see B12 for details; see also Girardi et al. 119971 Girardi 
et al. 2010). Moreover, we investigated the results of the DS-test 
obtained by changing the number of neighbors. 

The two panels of Fig.fTOlillustrate the significant results for 
the two kinematical indicators. This figure shows the spatial dis- 
tribution of the 87 cluster members. Each galaxy is indicated by 
a bubble, where the size of the bubble is related to the value of 
5, the deviation of the local kinematic estimator from the global 
cluster value. As for the 6y estimator, the substructure signif- 
icance increases up to > 99.9% c.l. when A n n =40. Figure [TUl 
(upper panel) shows that there are two regions of different lo- 
cal mean velocities at NE and SW, in agreement with the veloc- 
ity gradient discussed above, and likely corresponding to a low- 
and a high-velocity galaxy populations. When considering the 5 S 
estimator, the substructure significance is higher at smaller N nn . 
Figure [10] (lower panel) shows a NE region of low local velocity 
dispersions (N Dn =10, significance at the 96% c.l.). Probably, in 
that region, there is a minor mixing of the low- and high-velocity 
populations. 

We also applied the full 3D-KMM method. Considering 
GV1 and GV2 as a guess for an initial two-group partition, we 
find that the 3D galaxy distribution is well described by a parti- 
tion of 49 and 38 galaxies (KMM3D-1 and KMM3D-2 groups). 
The improvement over the sole 3D Gaussian is significant at the 
98% c.l. Figure [TT1 illustrates the distribution of KMM3D-1 and 
KMM3D-2 on the plane of the sky and the two Gaussians in the 
velocity distribution. Table [2] lists the main properties of these 
two groups. With this 3D test, the values of the velocity disper- 
sions of KMM3D-1 and KMM3D-2 are much higher than those 
of the corresponding groups obtained with the ID methods pre- 
viously described. 

We finally applied the Htree method (Serna & Gerbal [l996l) . 
This method performs a hierarchical clustering analysis and re- 
turns the relationship between cluster members based on their 
relative binding energies (see also, e.g., B12 and Girardi et al. 
1201 11 ). Figure [T_2] illustrates the results of the Htree method. In 
this dendogram the binding energy is reported in abscissa (in ar- 
bitrary units) and galaxy pairs and subgroups lie to the left (at 
lower energy levels). Going down from the top of the dendo- 
gram, we note the secondary system HT2 and the main system 
HT1. HT2 is a group at high velocity. HT1 is the main system 
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Fig. 10. Positions on the plane of the sky of 87 cluster members 
(marked by bubbles). Top panel: DS-test bubble plot for the de- 
viation dvj and A nn = 40 (see text for an explanation). Bottom 
panel: as above, but for the deviation 6 s j and A nn = 10. In both 
panels, thin/blue (thick/red) bubbles show regions with a local 
value lower (higher) than the global one. 



and shows a low-velocity substructure (HT12), while its core 
(HT11) is dominated by the BCG The spatial distributions of 
galaxies of HT1, HT2, and HT12 are plotted in Fig. [13] When 
using the results of the Htree method as seeds for the 3D-KMM 
method, we do not find any two or three-group partition signifi- 
cantly better than the sole 3D-Gaussian. 
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Fig. 11. 2D distribution of 87 cluster members. Galaxies with 
smaller (larger) symbols have higher (lower) radial velocities. 
Big blue circles and red crosses indicate galaxies of KMM3D-1 
and KMM3D-2. The insert plot shows the same velocity dis- 
tribution of Fig. [5] and the two Gaussians corresponding to 
the mean velocity and velocity dispersions of KMM3D-1 and 
KMM3D-2 (blue thin line and red thick line, respectively). 



4. X-ray morphological analysis 

We studied the morphological properties of the ICM by using 
archival X-ray data taken with Chandra ACIS-S (exposure ID 
#906, total exposure time 58 ks). We reduced the data using the 
package CIACQ (ver. 4.2) on the chip S3 in a standard way (see, 
e.g., Boschin et al. 120041 . 

The reduced image (see Fig. [T4l reveals that the ICM ex- 
hibits a regular morphology. From a quantitative point of view, 
this result is supported by the power-ratio (Buote & Tsai 1996) 
analysis performed by Hart (2008). Moreover, we used the task 
CIAO/Wavdetect on chip S3 to perform a wavelet multiscale 
analysis and found no evidence of multiple clumps in the ICM. 

We also used the CIAO package Sherpa to fit an elliptical 
/?-model profile (e.g. Boschin et al. [2004) to the X-ray pho- 
ton distribution (after removing of the pointlike sources found 
with CIAO/Wavdetect). The best-fit model is characterized by 
a centroid position located ~4" east of the BCG and a core ra- 
dius ro = 47.3±1.2 arcsec (221+6 hzl kpc). The best-fit val- 
ues for other parameters are: e = 0.213+0.007 (ellipticity), 
a = 1.74+0.04 (power law index), and 0=149.8+1.1 degrees 
(angle of ellipticity). 

The reduced CSTAT statistic (Cash [19791 of the fit is 1.22. 
Thus, the elliptical beta model is an adequate description to the 
data. However, we checked for possible deviations in the X- 
ray photon distribution from the above model by computing the 
model residuals. We find a deficiency of photons in a dumbbell- 
shaped region extending along the SSE-NNW direction in prox- 
imity of the X-ray centroid position (see Fig. [TBI. Regions with 



4 see http://asc.harvard.edu/ciao/ 




Fig. 12. Results obtained with the Htree method (see text). In this 
dendogram, the horizontal axis reports the binding energy while 
the vertical axis shows the IDms of the member galaxies (as in 
TableQ]). 

positive residuals elongated in the NE-SW direction are also 
found all around the cluster center. 



5. Discussion and conclusions 

Our estimate of the velocity dispersion (<x v = 1302^j 7 km s~ ] ) 
agrees with those of Patel et al. (2000) and Irgens et al. (2002 ) 
computed on a much smaller galaxy sample. This high value of 
cry predicts kTx = 10.3*|j keV (assuming energy equipartition 
between galaxies and gas energy per unit mass), and thus agrees 
with the high measured X-ray temperature kTx = 7-9 keV 
(from Chandra data, see Baldi et al. 2007 Bonamente et al. 2008 
and Ehlert & Ulmert 120091 . Both cr v and T x suggest that A1995 
is a massive galaxy cluster. 

In the hypothesis of dynamical equilibrium (but see in the 
following) and typical assumptions (cluster sphericity, galaxies 
and mass following identical distributions), we computed virial 
global quantities. By considering Girardi & Mezzetti (200T| see 
also B12 for details), we obtained a measurement of the mass 
within the virial radius R vil : M(< R vil = 2.7 /i^Mpc) = 3.0^ 
xlO 15 h'l M Q . 

5. 1. Cluster structure and mass 

Substructure is detected using several analyses of the cluster 
galaxy population. Our optical analyses indicate the presence of 
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Fig. 13. 2D distribution of 87 cluster members. Big blue circles 
and red crosses indicate galaxies of HT1 (main system) and HT2 
(high-velocity subcluster). The insert plot shows the same veloc- 
ity distribution of Fig. [5] (dashed line) and that of HT1 (blue thin 
line) and HT2 (red thick line). The subcluster HT12 of HT1 is 
indicated by small blue crosses in the main plot and this causes 
the low-velocity tail in the insert plot. 



two main subclusters aligned in the NE-SW direction causing 
the velocity gradient towards the NE direction and separated by 
~ 1.5', i.e. a projected linear distance D ~ 0.4 hz}. Mpc. 

The two subclusters, as determined through ID- or 3D- 
KMM methods, have comparable velocity dispersions within the 
errors, and in both cases they are likely two massive systems. 
The strong uncertainty in the subcluster membership reflects on 
the recomputation of the system mass by summing the masses 
of the two subclumps; i.e., M sys ~ 2-5 xlO 15 M Q , where the 
low (high) value is computed for the ID (3D) case with a rest- 
frame velocity separation of AVrf _los ~ 2000 km s _1 (~ 600 
km s _1 ; see Table|2|. Both ID and 3D methods have their draw- 
backs. As for the ID case, the two peaks are not clearly sep- 
arated in the velocity distribution. In the 3D case, our spectro- 
scopic sample is not spatially complete and less extended than 
the supposed /?2oo- We thus suggest that intermediate values are 
closer to the real one. This leads to a mass value not far from the 
virial mass previously computed considering the global velocity 
dispersion. 

Moreover, we also have indications of a more complex struc- 
ture: the 2D-DEDICA analysis detects a third, minor central 
galaxy peak aligned along the NE-SW direction. Finally, the 
Htree method does not detect a clear bimodal structure. 

For the X-ray data analysis, we reaffirm the existence of an 
isophotes elongation. More quantitatively, we find isophotes are 
elongated with a position angle PA = #-90°=59.8°(where PA is 
measured in a clockwise direction from the north; see Sect. 3J. 
This PA agrees with the direction of the velocity gradient and 
is slightly larger than the PA suggested by the direction defined 
by the two main optical subclumps. This scenario is somewhat 
suggestive of a cluster merger. 




Fig. 14. Chandra image of A 1995 (smoothed X-ray emission in 
the energy band 0.3-7 keV, direct sky view with north up). The 
field of view is 8.5'x8.5'. 



Moreover, while previous studies did not find any direct ev- 
idence of substructure in the X-ray emission (Ota & Mitsuda 
2004; Pedersen & Dahle 2007), our detailed analysis of Chandra 
data suggests the presence of multiple clumps. The elliptical (3- 
model residuals resemble the residuals computed by Girardi et 
al. (|20l0) from the fit of the X-ray photon distribution of the 
cluster Abell 2294 (see their Fig. 12) with a simple /3-model. 
For Abell 2294, Girardi et al. (2010) explain the residual image 
by the presence of two very closely projected subsystems. Our 
Fig-EJsuggests a similar scenario for A1995, with an excess of 
X-ray emission (positive residuals) in the NE-SW direction. The 
main central X-ray peak is close to the BCG This peak and a 
SW one are offset with respect to the two main optical peaks and 
located between them, thus resembling the usual situation of a 
bimodal merging just after the core-core passage with the colli- 
sional X-ray components slowed down with respect to the col- 
lisionless galaxy clumps. The NE X-ray peak lies between the 
main optical peak and the third (minor) optical peak, thus sug- 
gesting that a third subcluster intervenes in the merger. A fourth 
X-ray peak (excess of positive residuals) at the NW is less sig- 
nificant and is not considered. Summarizing, our results confirm 
the active dynamical state of A 1995 from the thermodynamic 
point of view. 

5.2. Merger kinematics and diffuse radio sources 

Considering only the two main subclusters, the merger can be 
followed through the analytic method of the bimodal model 
(Beers et al. 119821 see also Boschin et al. 2010). For the val- 
ues of the relevant parameters (D; M sys ; AVrf^osX we use d the 
values reported above taking both the cases corresponding to 
ID- and 3D-KMM analysis into account (hereafter case- ID and 
case-3D). We assumed t = 0.2 Gyr for the time elapsed from 
the core crossing (a typical value for clusters hosting radio ha- 
los; e.g. Markevitch et al. 120021 Girardi et al. l2008l . Figure [TBI 
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Fig. 15. The smoothed X-ray emission of A1995 (direct sky 
view with north up). Black (white) contours show the positive 
(negative) elliptical yS-mode residuals. Big crosses indicate the 
two main optical clumps detected with the 2D-DEDICA method 
(2D-NE and 2D-SW in Table©. Small cross indicates the third 
densest peak we detect (2D-NENE). 



illustrates the solutions of the bimodal model, where the mass 
of the system M sys is plotted against a, the projection angle be- 
tween the direction defined by the subclumps and the plane of 
the sky. In both case- ID and case-3D we find a bound outgoing 
solution (BO, see Fig. [Toll compatible with a ~ 50°. This value 
agrees with the strong evidence of substructure in the 3D analy- 
ses, since 3D substructure is more easily detected at intermediate 
angles (Pinkney et al. 119961 . Note that the assumption of signif- 
icantly longer times t would lead to larger, unlikely, angles of 
view (e.g. > 70° for t > 1 Gyr). 

As a comparison with the results in the radio band, 
Giovannini et al. (|2P09) show that, as usual, the radio halo of 
A 1995 and the X-ray emitting ICM occupy the cluster volume in 
the same manner (see their Fig. 8 and our Fig.[TJ. We find that the 
elongation of the radio halo roughly agrees with the direction of 
the merger, as detected in several clusters, e.g. the Bullet Cluster 
(1E0657-56; Markevitch et al. 120021 and Ab ell 520 (Girardi et 
al. 120081 : but see Abell 523 (Giovannini et al. [20TTT >. 

In conclusion, A1995, for its high mass and cluster merg- 
ing evidence, is not atypical among clusters hosting radio halos. 
The remaining puzzling point concerns the dark matter distri- 
bution, which is quite circularly symmetric with respect to the 
galaxy distribution (Dahle et al. 120021 Holhjem et al. l2009l l. the 
latter indicating a NE-SW preferential direction. Owing to the 
assumed collisionless nature of both galaxies and dark matter 
particles, this is surprising evidence. Moreover, this disagrees 
with the results found for other extensively studied clusters, such 
as the Bullet Cluster (e.g., Markevitch et al. 120021 or clusters 
CL 0152-1357 and MS 1054-0321 (Jee et al. I2005al l2005bl . 
where a coincidence is found between the galaxy and dark mat- 
ter distributions. Nevertheless, the peculiarity of A1995 is not 
so extreme as in other clusters, such as CL 0024+17, where ev- 
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Fig. 16. System mass - projection angle diagram of the ana- 
lytic bimodal model applied to the two subclusters. Thin/red 
(thick/blue) lines refer to the case-ID (case-3D; see text). In par- 
ticular, bound (solid curves) and incoming (collapsing; BI a and 
Bib), bound outgoing (expanding; BO), and unbound (dashed 
curves) outgoing (UO) solutions are shown. The horizontal lines 
give the computed values of the system mass for the two cases. 
The dotted curves separate bound and unbound regions (above 
and below the dotted curves, respectively). 

idence was found for a dark matter ring (Jee et al. 120071 1 with- 
out a obvious galaxy counterpart (Qin et al. 120081 but see also 
Ponente & Diego 1201 11 who suggest that dark matter ringlike 
structures could be due to systematics in lensing reconstruction). 
In any case, the observed dichotomy between the dark matter and 
galaxy distributions makes A 1995 an appealing target for future 
studies. 
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